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Preface 


The  purpose  of  this  thesis  was  to  study  the  propagation 
of  magnetostatic  waves  on  a  YIG  ring.  This  topic  was  chosen 
because  of  my  interests  in  experimentally  oriented  research 
and  magnetostatics.  I  saw  the  opportunity  to  expand  on 
someone's  initial  effort,  with  ideas  and  techniques  of  my 
own.  This  thesis  topic  provided  me  the  chance  to  actually 
build  electronic  devices  and  test  them.  I  believe  the 
following  findings  and  conclusions  will  be  useful  for  those 
participating  in  research  in  the  magnetostatic  wave  area. 

I  am  most  deeply  grateful  to  Mr.  James  Sethares  of  the 
Rome  Air  Development  Center's  Electromagnetic  Sciences 
Division  for  sponsoring  the  thesis,  in  providing  samples  for 
testing,  advice  during  experimentation,  and  mostly  his 
encouragement  and  guidance  during  my  first  years  as  a 
scientist  and  engineer. 

A  word  of  thanks  also  goes  to  Marc  Calcatera  of  the 
Avionics  Laboratory  at  Wright-Patterson  AFB  for  his  help 
during  the  early  experiments.  Finally,  I  am  extremely 
thankful  for  the  continuing  patience  and  assistance  of  my 
faculty  advisor,  Capt.  Roger  Colvin. 


Robert  E.  Floyd 
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Abstract 


A  unique  test  apparatus  was  designed  in  order  to  inves' 
tigate  the  effect  of  launching  magnetostatic  surface  waves  at 
different  angles  on  a  yttrium-iron-garnet  (YIG)  single 
crystal  film.  Many  characteristics  of  the  magnetostatic 
waves  were  measured,  including  amplitude,  time  delay,  and 
beam  steering.  There  was  slight  dependence  of  these 
characteristics  on  the  launch  angle.  A  second,  more 
sensitive  technique  using  the  YIG  films  as  delay  lines  in  an 
oscillating  circuit,  was  found  to  verify  the  influence  of 
magnetocrystalline  anisotropy  on  the  magnetostatic  surface 
wave.  A  calculated  shift  in  frequency  of  the  oscillator  was 
found  to  be  in  agreement  with  previous  calculations  and 
measurements  of  anisotropic  constants. 

In  addition  to  these  findings,  an  investigation  of  the 
influence  of  a  metallization  spaced  closely  to  the  propa¬ 
gation  path  was  performed  due  to  its  influence  on  the  above 
measurements.  This  metallization  altered  the  velocity  of  the 
magnetostatic  wave  and  demonstrated,  for  the  first  time,  the 
refraction  of  a  magnetostat i c  wave.  This  technique  may  prove 
useful  for  instantaneous  frequency  selection  in  electronic 
countermeasure  equipment. 


I .  Introduction 


Background 

The  work  of  Damon  and  Esbach  in  the  late  1950 's  and 
early  1960's  provides  the  theory  on  which  much  of  magneto¬ 
static  wave  technology  is  based  (Damon,  I960).  Their 
derivations  show  the  existence  of  magnetostatic  surface  wave 
modes  in  ferromagnetic  materials.  This  behavior  is  similar 
to  spin-wave  propagation  in  an  infinite  media  when  no 
boundaries  exist.  With  the  introduction  of  boundaries,  and 
as  the  excitation  frequency  is  increased,  the  spin-wave 
behavior  changes  from  one  that  is  distributed  throughout  the 
volume  of  the  medium,  to  one  that  is  confined  to  the  surface. 
Analysis  indicates  that  these  two  types  of  modes  are  unique, 
frequency-dependent  solutions  to  Maxwell's  equations,  and  the 
equation  of  motion  of  the  slab's  magnetization  (Damon, 
1960:309).  Surface  wave  propagation  is  initiated  by  external 
RF  magnetic  fields.  The  energy  of  the  surface  mode  is 
concentrated  along  the  surface  where  propagation  takes  place, 
with  the  propagation  direction  orthogonal  to  an  applied  DC 
field. 

In  1968.  experiments  by  Brundie  and  Freeman  verified 
some  of  Damon  and  Esbach's  theory.  Brundie  and  Freeman  were 
able  to  excite  the  surface  wave  mode  on  a  thin,  single 
crystal,  yttrium-iron-garnet  (YIG)  slab  where  the  magneto¬ 
static  surface  wavelengths  were  much  smaller  than  the  slab 
dimensions  (Brundie,  1968).  Results  closely  agreed  with 


Damon  and  Esbach's  theory  in  terms  of  time  delay,  propagation 
direction,  and  dispersion.  The  surface  wave  was  observed  to 
be  dispersive,  and  the  wave  velocity  decreased  as  the 
excitation  frequency  was  increased  (Brundle,  1968:427). 

Many  signal  processing  applications  have  been  developed 
using  magnetostatic  surface  wave  (MSSW)  technology.  Tunable 
microwave  oscillators  (Miller,  1976)  and  variable  delay  lines 
operating  in  the  microwave  spectrum  have  been  developed. 

Many  are  based  upon  the  theory  of  Damon  of  Esbaoh,  the 
experiments  of  Brundle  and  Freeman,  and  others  (Collins, 

1977).  This  list  of  devices  should  also  include  the  proposal 
of  Newburgh,  Blacksmith,  Budreau,  and  Sethares  of  the  Air 
Force  Cambridge  Reserch  Laboratory  (Newburgh,  1974).  These 
researchers  proposed  the  development  of  a  small,  highly 
sensitive  and  lightweight,  non-inertial  rotation  rate  sensor 
using  MSSW  propagation.  Using  an  interferometric  sensing 
scheme,  the  rotation  rate  would  be  detected  by  measuring  the 
phase  shifts  of  magnetostatic  surface  waves  on  a  rotating  YIG 
ring.  In  1974,  Sethares  showed  that  MSSW  propagation  indeed 
could  be  achieved  on  a  YIG  thin  film  ring  grown  on  a 
gadolinium-gallium-garnet  (GGG)  substrate  (Sethares,  1975) 

(see  Figure  1.1).  The  use  of  the  thin  film  ring  as  a 
rotation  sensor  proved  questionable  because  the  experimental 
results  exhibited  large  variations  in  coupling  efficiency  as 
the  ring  was  rotated  about  its  axis.  Attributing  the  results 
to  imperfections  of  the  growth  in  the  thin  YIG  film,  Sethares 


suggested  the  use  of  a  thick  film  ring.  This  ring  would  be 
machined  from  a  single  crystal  YIG  boule,  so  that  the  crys¬ 
talline  imperfections,  which  occur  during  the  liquid  phase 
epitaxial  growth  process  of  a  thin  film,  could  be  eliminated 


Figure  1.1  YIG  Film  on  the  Outer  Surface  of  a  GGG  Ring 
Substrate  (Sethares.  1975) 

An  evaiuati'  n  of  the  thick  film  ring  geometry  was  done 
by  Poturalski  in  a  master's  thesis  in  1979  (Poturalski, 

1979).  Large  coupling  variations  were  significantly  reduced 
and  Poturalski  was  able  to  propagate  surface  waves  around  th* 
ring,  but  was  unable  to  achieve  satisfactory  rotation  rate 
sensing  sensitivity.  Some  unexplained  variations  were  still 
present . 

Problem  Statement  and  Scope 

Poturalski  observed  a  non-uniformity  of  coupling  and 
dispersion  as  a  function  of  lattice  orientation  at  the 


surface  wave  launch  site  ( Poturalski ,  1979:80).  He  suggested 
that  an  analysis  be  performed  on  the  effect  of  the  orien¬ 
tation  of  the  crystal  lattice  at  the  surface  wave  launch  site 
and  the  propagation  direction.  This  suggestion  serves  as  the 
premise  for  the  present  investigation. 

Since  the  rotation  rate  systems  proposed  by  Sethares 
use  the  characteristics  of  magnetostatic  surface  wave  (MSSW ) 
propagation  to  sense  rates,  the  experimental  analysis  per¬ 
formed  on  the  1-inch  YIG  samples  in  this  thesis  pertains  to 
the  characteristics  of  surface  modes  only.  To  simplify  the 
geometry,  the  orientation  effects  of  crystalline  direction 
were  studied  on  flat,  thin  films  of  YIG  deposited  on  a  GGG 
surface . 

Several  characteristics  of  the  surface  mode  were 
investigated.  Insertion  loss  measurements  between  the 
transmitting  and  receiving  transducers  were  made  to  observe 
the  effects  of  signal  power.  Time  delays  of  pulsed  RF 
signals  were  included  because  of  the  MSSW"  s  technological 
importance  to  variable  time  delay  devices.  All  character¬ 
istics  were  measured  as  a  function  of  a  rotation  angle  (HI 
which  is  measured  relative  to  the  crystal  lattice  orien¬ 
tation.  In  addition  to  varying  the  rotation  angle,  the 
receiver  transducer  was  varied  in  both  axes  of  the  YIG  film 
plane  as  shown  in  Figures  1.2  and  3-3  in  order  to  obtain  pro¬ 
file  measurements  of  the  wave  as  it  travelled  across  the 
film.  The  ability  to  independently  vary  these  parameters  is 


facilitated  by  the  coupling  probes  designed  and  developed  by 
Poturalski  for  the  rotating  YIG  ring  (Poturalski,  1979:39) 
(see  Figure  3.2). 


i; 


Figure  1.2  Top  View  of  Apparatus  Used  in  the  Investigation 

11  though  the  Poturalski  probes  were  indispensable  to 
the  rotation  measurements,  they  were  suspect  for  the  intro¬ 
duction  of  measurement  errors.  An  evaluation  of  the  probes 
detected  a  steering  effect.  The  metal  cover,  which  was  used 
to  prevent  electromagnetic  radiation  feedthrough,  was  thought 
to  be  the  cause  of  this  phenomena.  An  experiment  to  confirm 
the  influence  of  a  metallization  spaced  closely  to  the 
propagation  path  was  performed.  These  measurements  produced 
the  most  rewarding  findings  of  this  report. 


Overview 


For  the  reader  who  is  unfamiliar  with  MSSW  technology, 
the  derivation  of  the  electromagnetic  field  equations  for  the 
magnetostatic  modes  and  surface  wave  propagation  is  discussed 
in  the  second  chapter.  This  will  enable  the  reader  to  under¬ 
stand  the  interactions  of  externally  applied  magnetic  fields 
and  ferromagnetic  materials.  Also,  a  pertinent  discussion  of 
crystallography  and  its  importance  to  the  investigation  is 
included . 

The  design  and  development  of  experimental  equipment 
and  procedures  is  described  in  the  third  chapter.  This 
chapter  discusses  the  ferromagnetic  samples,  the  probes  used 
for  transducers,  the  experimental  test  fixture,  and  two  test 
configurations  used  for  gathering  magnetosta ; i c  character¬ 
istics  versus  various  surface  wave  launch  angles. 

Analysis  of  the  results  obtained  from  the  experiments 
are  discussed  in  the  fourth  and  fifth  chapters .  Analysis  of 
radio  frequency  (RF)  pulse  measurements  revealed  a  slight 
dependence  of  these  magnetostatic  characteristics  on  the 
surface  wave  launch  angle.  A  second,  more  sensitive  tech¬ 
nique  using  an  oscillating  circuit,  was  found  to  verify  the 
influence  of  magnetocrystalline  anisotropy  on  the  magneto¬ 
static  surface  wave.  A  calculated  shift  in  frequency  of  the 
oscillator  was  found  to  be  in  agreement  with  previous  calcu¬ 
lations  and  measurements  of  anisotropic  constants. 

It  was  found  during  the  initial  experiments  that  the 
probes  had  an  influence  on  the  pulse  and  oscillator 
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measurements.  In  the  fifth  chapter,  an  investigation  of  a  . 

t 

metallization  spaced  closely  to  the  propagation  path  was  j 

i 

performed  to  determine  its  influence  on  the  magnetostatic  ! 

i 

wave.  This  metallization  altered  the  velocity  of  the  ! 

I 

magnetostatic  wave  and  demonstrated,  for  the  first  time,  the  j 

refraction  of  a  magnetostatic  wave. 

The  refraction  technique  may  prove  useful  for  i 

instantaneous  frequency  selection  in  electronic  counter-  j 

measure  equipment.  This,  and  other  recommendations  for 
further  study,  along  with  concluding  remarks,  are  found  in 
the  last  chapter. 
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II .  Magnetostatic  Surface  Wave  Theory 


Background 

A  description  of  the  interactions  between  internal 
magnetic  moments  and  externally  applied  magnetic  fields  will 
best  explain  the  magnetostatic  behavior  in  ferromagnetic 
materials,  such  as  yttrium-iron-garnet  (YIG). 

Sparks  classically  explains  magnetostatic  behavior  in 
terms  of  magnetic  moments  due  to  electron  spins  within  a 
crystal  lattice  (Sparks,  1964:5).  In  ferromagnetic  mate¬ 
rials,  neighboring  spins  and  their  associated  magnetic 
moments  tend  to  be  aligned  parallel  to  one  another.  By 
integrating  over  a  volume  of  these  magnetic  moments,  one 
obtains  a  total  magnetic  moment  which  can  be  perturbed  by 
external  fields. 

Next,  considering  the  spin  system  when  a  DC  field  is 
applied,  a  change  is  evident.  Intuitively,  this  field  has  an 
effect  on  an  isolated  one-spin  system,  and  thus,  on  the  total 
magnetic  moment.  Further,  one  can  show  how  this  spin  system 
is  related  to  the  magnetic  moment.  A  torque  is  produced  in 
the  presence  of  the  DC  field  which  causes  a  precession  of  the 
magnetic  moment  around  its  unperturbed  position.  The  rate  of 
this  precession  is  referred  to  as  the  Larmor  frequency,  and 
the  magnetic  moment  rotates  at  some  precession  angle  that  is 
measured  with  respect  to  the  unperturbed  magnetic  moment 
direction  (Figure  2.1). 


DC  FIELD  =  H 


*  -  -  —  ■ 
Figure  2.1  One  Spin  System  (Sparks.  196^) 

Sparks  discusses  the  magnetic  moment's  reaction  to  an 
RF  field  applied  orthogonally  to  the  DC  field.  When  the 
frequency  of  the  RF  field  is  set  to  the  precession al  rate  of 
the  spin  system,  a  two-fold  change  is  seen  in  the  magnetic 
moment  vector.  The  magnitude  of  the  magnetic  moment  vector 
and  the  precession  angle  are  both  increased  (Sparks,  196^:6) 
These  changes  are  represented  by  Mi  and  M2  in  Figure  2.1. 

For  any  rotating  body ,  the  rate  of  change  of  the 
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angular  momentum  J  equals  the  applied  torque  T: 


The  precession  of  rotation  is  at  a  rate  sufficient  to 


conserve  the  initial  angular  momentum.  Another  relationship 
of  the  applied  torque  involves  the  magnetic  moment: 

(2> 

where  %  is  the  gyromagnetic  ratio. 

The  torque  resulting  from  subjecting  a  magnetic  moment 
m  to  a  magnetic  field  H  is  given  by: 

T=5xH,  (3) 


Combining  (1)  and  (3)  with  (2)  yields 

=  %  (  **  X  H)  (jj) 

at 

which  is  the  equation  of  motion  for  the  magnetization 
(Emptage,  1982:266). 

It  is  noted  that  a  demagnetization  appears  when  this 
single  system  is  expanded  to  include  more  spins.  This 
opposition  is  caused  by  a  reactionary  field  which  acts  to 
demagnetize  the  applied  DC  field.  Two  properties  which 
linearly  relate  to  this  demagnetization.  H^,  are  a  demagnet¬ 
ization  factor,  Nj,  which  is  a  function  of  the  material’s 
geometry  and  the  saturation  magnetization,  MnMs,  of  the 
ferromagnetic  material  (Soohoo,  1960:55). 


This  alters  the  net  effective  field  which  now  becomes 

H*  =  -  Hj  (6) 

where  Ha  =  applied  field. 

The  preceeding  discussion  does  not  include  all  physical 
phenomena  that  will  affect  magnetostatic  behavior.  Limita¬ 
tions  in  the  model  include  omission  of  temperature  effects, 
impurities,  and  surface  defects.  These  issues  are  important 
and  several  will  be  discussed  later,  but  despite  these 
effects,  the  equations  provide  a  fundamental  basis  for 
discussion  of  the  magnetostatic  surface  wave  mode  of 
propagation . 

Magnetostatic  Surface  Waves 

Now,  with  a  basic  understanding  of  the  internal 
magnetic  characteristics  of  a  ferrite  such  as  YIG,  we  can 
place  boundary  conditions  on  the  ferrite  and  introduce  a 
current  transducer  which  provides  an  external  magnetic  field 
signal  we  would  like  to  propagate. 

This  phenomena  is  best  analyzed  by  s imu 1 taneous 1 y 
solving  Maxwell's  equations  and  the  gyromagnetic  equation, 
along  with  the  appropriate  electromagnet i c  boundary 
conditions.  Figure  2.2  illustrates  the  position  of  the 
boundaries.  The  conductor  between  Regions  I  and  II 
represents  the  transducer.  Region  III  is  the  YIG  film. 

Regions  I,  II,  and  IV  are  air  surrounding  the  film. 
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of  the 

MSSW  System 

l's  equations  are  satisfied  in  each  of  the 
four  regions: 
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For  MSSW  propagation  in  the  y  direction. 


HT  =  E„'Et'  o  03) 

and  using  the  magnetostatic  approximation  (Morgenthaler , 
1983:135),  we  assume  the  time  rate  of  change  of  D  is  approxi 
mately  zero  such  that. 


_  uo£  t  ¥  &  °  (1^) 

H  ' 

Maxwell's  equations  for  our  magnetostatic  system  are  now: 
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and 
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The  next  step  is  to  satisfy  the  magnetization  motion 


equation:  a  linear  form  is  considered,  such  that 


B 


» 


For  Region  III,  the  film's  anisotropy  requires 


(  18) 


(  19) 


whereas  for  Regions  I,  II,  and  IV,  we  have  the  unit  matrix. 


Using  Emptage's  notation  (Emptage,  1982:  5122),  the  relative 
permeabilities  are: 


(20) 

=  SL.  (  -TL  ) 

with 

(21  ) 

-A-  '  ‘~l  /  4  tr  )i 

(22) 

a-m  -  ney4irms 

(23) 

Now.  since  Dx  =  Dy  =  0,  and  there  is  no  z  dependence 

Equation  (10)  is  satisfied,  and  it  is  noted  that: 

of  D. 

D  =  &X-E 

so 

(24) 

D_  r  E2 

(25) 

with  t'p  =  1  for  Region  I,  II,  and  IV:  and  cp  -  13.0  (Lax, 
1962:703)  for  Region  III. 

Under  the  magnetostatic  approximation,  the  H  field  may¬ 
be  defined  in  terms  of  a  magnetic  potential  function  where 


In  order  to  solve  equations  (15-18)  with  such  a  potential,  we 
assume  a  solution  of  the  form 


(27) 


with  j  =  I,  II,  III,  IV  denoting  the  region  oi  interest.  It 
will  be  shown  that  aj  is  positive  and  is  determined  by 
equations  (15-18).  The  constants  Aj  and  Bj  are  determined  by 
the  boundary  conditions. 

By  assuming  a  time-invariant  solution  of  equation  (27). 
the  variable  t  can  be  suppressed.  Then,  by  inserting  the 
assumed  solution  into  equation  (26),  we  get 


in  all 
(29). 


H3j  -  -  k  (A, 

regions  with  j  =  I, 
From  equat ion  (  1 8 ) 

Bxj  ~j*c  ^  *j 


(28) 

(29) 

II,  III.  IV  for  equations 

(28)  and 

(30) 

and  for  the  YIG  region, 

TIL 


b»_-aA  e'^)e-ik* 


V  *  V  * 


13a  Snr£  )e  * 

^  ^••K^V^e-’r02) 

In  order  to  determine  aj,  we  look  for  solutions  to 
equation  (17).  From  equation  (30)  we  can  get  the  following 
for  the  non-YIG  regions, 

=  A,  «-j  ( Aj  e,jX+  t>;  *'**)&'  t8'  j  >1,1,11  <33) 

ix  7 

23»L  id  (Aj  «?*  4  -B,-  e  *j *) e-! *>  ” » 

H 

noting  that  if  aj  =  |K  ,  then  equation  (17)  is  satisfied.  In 


Region  III 


2*54  .^,4  (A«e*“' 

\T2 

_ !E=  7VV« Bme'AmXX-ik^ 

^  -/A  (Ve‘m"‘  +'Bmea,,'')£-'k5!-  . 

Now  equation  (17)  will  be  satisfied  for  Region  III  if: 


(35) 


( 

/naI 

Where  PJ=R22/R11  so  that  am  =  pK,  in  order  to  satisfy 
equation  (17). 


( 37 ) 


Vk  -.VAv'l  '  \  .  V-  1*.)-*  fc.«  .A  iV '*■-**.  ~ .  •' .  *' *  •*".  **. *'■  **•  ■' 
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Using  equation  (16)  with  the  help  of  equation  (30),  we 


integrate  and  obtain. 


(38) 

)e-\ 

(39) 

Again,  note  everything  is  satisfied  if  am  =  p!K  ,  and 


B^e**Vk* 


(40) 


^  4  Bor a“m')e'ik*  («i ) 

+  e“" *-Bm e'"“*)  r k»  . 

Equations  (40)  and  (41)  are  equal  if  equation  (37)  holds 
true.  Therefore, 


ctj  =  1  (42) 

where  j  =  I ,  II,  IV,  and 

|  k|  (i4?) 

Having  aj  determined,  we  can  now  determine  Aj  and  Bj  via  the 
boundary  conditions. 

From  Figure  2.2.  it  can  be  observed  the  YIG  sample  is 
in  free  space,  and  there  are  no  ground  planes  to  consider. 

The  conductor  at  x=-g  is  where  the  current  is  introduced. 

The  H  field  will  not  be  continuous  across  this  boundary.  The 
remainder  of  the  boundary  conditions  are: 


2. 

@ 

x  =  -g 

Bx  is  continuous 

(45) 

3. 

§ 

o 

ii 

X 

Hg,  Bx  are  continuous 

(46) 

4. 

e 

X  =  t 

H,  Bx  are  continuous 

(47) 

5. 

§ 

x  :  +0C 

CD 

X 

ii 

o 

(43) 

From  V  B 

=  o, 

,  we  get  equation 

(17),  and  using  equations 

(26) 

and  (30) 

we 

obtain : 

/,^  +/«. 
/  dK  / 

,  o 
^”4 

(49) 

Thus,  equation  (49)  reduces  to  the  Laplacian,  when 
411=422’  which  holds  true  for  surface  waves.  This  informs  us 
that  there  is  one,  and  only  one,  solution  for  -v  with  the 
above  boundary  conditions  imposed  (Adler,  1969:521).  Using 
our  assumed  solution,  we  need  to  determine  A  i  and  B  j  .  jn 
order  for  our  solution  to  satisfy  all  K ,  we  integrate  (2 9 ' 
and  (30,1  over  K. 

Hjj  =  j  -:tLKe'K,“  ♦*, j*I,l,E 

-«o  j  \  O'j  ) 

"  j  "llc  (^nre^'K>X  !  =  !) 

~joo  '  "  ' 


(53) 


By  writing, 


°( i  —y^-xi.  ^ +  sy*  u. 


(54) 


;  s=  */\t\  (55) 

we  can  write  equation  (53)  as: 


We  can  now  apply  boundary  conditions  at  x  =  -  *,  where  we  can 
see  that  Bxj  =  0,  therefore,  in  Region  I, 


(57) 


and  at  x  =  -g,  BXI  =  BxIIf  and 


,f\T 

a-|U3-Bi[e,^3 

(58) 

at  x  =  0,  Bxii  =  BxIII  ,  Hyii 

:  =  Hy III,  SO 

~  Bj-  =  <3, ,  “ 

(59) 

Ar  +  \  =  + 

(60) 

at  x  =  t,  Bxji  =  B  x I v »  Hy ill  = 

Hy IV’  SO 

^  lk.lt.  _  ,  ||clt.  -r\  -/klt 

(61 ) 

1  /»lt  r,  -IKlt 

Vs  -  Bjz  e 

(62) 

at  x  =  -  xBxIV  =  and 

M 

0 

• 

(63) 
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From  equations  (57)  and  (58) 

V-  °  • 

From  equations  (61),  (62),  and  (63) 

-  ° 


Pjj.  =Ane°VVKlt+  ■gare<|-'s)|llt 


From  equations  (59)  and  (60) 


*  '4  [  Affl  ( 1  +  af’) 1  sm  ( i  -<o] 

B*  =  &  [Am  0 -«.■)■<■  4  C I  +  o] 

Letting 

V  *  (fLiU 

'  -  I  / 

we  can  write  equation  (61)  using  equations  (66)  and  (67) 

Next,  by  defining 

Y  =  + 

equations  (68)  and  (69)  become 


Using  equation  (65)  the  results  can  be  expressed  as 


4  =  *3  (X-Y^) 


Now  all  constraints  are  defined  in  terms  of  Am- 

One  more  boundary  condition  is  needed  to  solve  for  Am. 

This  last  condition  introduces  the  current  x  =  g.  That  is, 


(78) 

The  J(y)  term  is  a  prescribed  surface  current  distribution 
function.  In  order  to  account  for  the  generation  of  all 
possibilities,  equation  (78)  is  written  as  an  integral  over 
K.  That  is, 

=  j.1,1  (79) 


at  x  =  g,  and 


''5  KfAjr0e'“5t 

-A^V't  -  Sr«)e  "‘ls)e"k*Jk 


(80) 


•  1  * 

Multiplying  both  sides  by  ei*  y  and  integrating  with 
to  y  from  -«  to  <*,  the  result  is 

'Ii^e"k3d 

_>o©  ~0C 

By  noting  the  following  two  definitions,  the  Fourier 
transform  of  J(y)  is: 

c'k^<j^  -  «JV*0 

and  with  the  Dirac  delta  function: 


equation  (81)  reduces  to 

IOO 


I  <JT(0-  J  (y  k' 


Upon  evaluation,  the  result  is 


OV^-SHe'1'3 


respect 


(81) 


(82) 


(83) 


(814) 


(85) 


by  the  sifting  property  of  the  delta  function  (Marsden, 
1981:392). 


Using  the  constraints  found  from  the  boundary 


conditions,  equation  (85)  reduces  to 


Ztt  Z  / 

*(4je'luV''3 


and  upon  simplification,  bee  es 


&  Ap  k. e/'"^ j^|X-7e1"‘,JJ -X)e'tlJ-Ye'""s 


Therefore,  with 


Z£|!lX-Ye,,,i-X)e,M,-Te"MiI 


then 


TT  K  e 


^  *  z 


Using  substitution  and  equation  (7*0 

A  =  1 .  2L. 

^  2^  Z 

and  from  equations  (75),  (76),  and(77) 

B_  -  i  J*L.  X- 

1  InK  Z 

-  icJVl  (X-Ye?'"* 


Using  equation  (71),  and  upon  simplification 


P 


*  W  J1  ^  .  Sl 

IT  K.  Zi 


IHt 


(9* 


From  equations  (65)  and  (67)*  these  constraints  can  now  be 
defined  as 

A  ~ 

( 9. 


A  =  ° 

*  '  JSL 


( I  ♦  W) 


«JJV> 


(9f 


tr  tc 


z 


From  equations  (50)  —  (53)  *  the  field  equations  can  be 
expressed  in  terms  of  the  current  density.  They  are 


-  I  -  fett* 


ik(>  _;k* 

e  «g.  *  ak 
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(91 


(9E 


u  rjw 

to  z^J„  ^ 
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•At  N 


IkI  JnqO+W 


where : 


V-  (°t+ 

X  =(><*,  )e^HM*  (/-^)e/lkU 

^[((X-T-'kl3)-X)e,Kl3-Ye,M( 


(  105) 


(  106) 


(  107) 


(  108) 


Thus,  all  of  the  constants  have  been  solved  for  in 
terms  of  the  current  density.  The  integral  field  equations 
summarized  in  equations  (97)-(104)  have  been  evaluated  by  the 
method  of  residues  (Sethares,  1976:14).  These  equations  are 
consistent  with  those  of  Sethares  (Sethares.  1976),  and  may 
be  used  as  a  starting  point  for  calculating  the  power  flow 
for  MSSW’s,  which  may  be  useful  in  future  investigations. 

From  Figure  2.2,  it  is  observed  that  the  magnetostatic 
wave  propagates  in  the  y-direction  with  an  amplitude 
dependence  on  x  (Sethares,  1975:13).  The  DC  magnetic  biasing 
field  is  in  the  z-direction. 


iV 


*.Y?  » 


I 

v  Since  the  wave  is  guided  along  the  surface  perpen¬ 

dicular  to  the  magnetic  biasing  field,  it  is  given  the  name 
surface  wave.  Experiments  by  Brundle  and  Freedman  verify  the 
direction  of  propagation  and  conclude  that  a  surface  wave  is 
indeed  present  (Brundle,  1968:1090). 

Another  characteristic  of  a  MSSW  is  that  it  is 
dispersive.  This  means  that  the  velocity  of  the  wave  is  a 
function  of  frequency.  This  relationship  is  useful  for 
interpreting  some  of  the  results  from  the  time  delay 
measurements ,  and  it  is  shown  in  the  followirg  equation  and 
Figure  2.3  (Sethares,  1975:14)  that: 

(f)  '  He  +  (^rVfCl  •  <'09) 

I 

where  k  =  MSSW  wave  number 
t  =  thickness  of  YIG 
w  =  MSSW  angular  frequency 
^  =  2n  (2.8xlo6)Hz  Oe  =  gyromagnetic  ratio. 


Crystal lography 

In  his  thesis.  Poturalski  (Poturalski,  1979:31),  recom¬ 
mended  the  investigation  of  crystalline  structure  effects  on 
MSSW  propagation.  What  follows  is  a  brief  description  of 
basic  crystallography  and,  specifically,  the  crystalline 
structure  of  yttrium-iron-garnet,  which  is  a  popular  propa¬ 
gation  media  for  MSSW.  This  analysis  provides  information  on 
how  the  orientation  angle  influences  the  magnetostatic  wave. 


'JL  +  JL 


Vo 
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Figure  2  3  Dispersion  Relation  for  MSSW  (Sethares,  1975:18) 

The  various  planes  that  pass  through  a  crystal  may  be 
described  in  terms  of  the  reciprocal  of  the  intercepts  of 
that  plane  with  the  crystallographic  axes.  These  are  writ’. e 
as  (h  k  1),  with  h.  k,  1  usually  expressed  as  the  smallest 
possible  set  of  integers  having  the  same  ratio.  This 
notation  format  is  known  as  the  Miller  index  (Wood.  19t?:8). 
A  bar  over  an  index  number  denotes  a  negative  intercept. 

When  one  refers  to  a  particular  crystalline  direction, 
the  format  is  to  enclose  the  Miller  index  in  brackets.  For 
example,  the  notation  Mill  is  the  direction  of  a  vector' 
normal  to  the  surface  of  a  (111)  plane. 

It  is  now  of  interest  to  introduce  the  structure  of 
YIG.  The  samples  to  be  investigated  are  all  thin  films  growr 
on  a  gadolinium-gallium-garnet  substrate.  The  films  were 


grown  by  liquid  phase  epitaxy;  thus,  the  YIG  will  take  the 
same  single  crystal  structure  as  the  GGG.  All  samples  tested 
had  (111)  surfaces  ((1111  direction  is  normal  to  the  surface 
of  the  sample). 

By  looking  at  this  structure  from  a  top-down  view,  as 
in  Figure  2.6,  one  can  see  that  the  lattice  structure 
exhibits  three  high  symmetry  directions. 
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To  test  for  crystal  orientation  dependence,  it  was 
necessary  to  construct  a  suitable  test  fixture  to  rotate  the 
samples.  This  experimental  apparatus  is  discussed  in  Chapter 
Three,  along  with  the  test  circuit  configurations. 
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Ill .  Experimental  Apparatus  and  Procedures 

Sample  and  Probes 

The  films  used  in  Sethares'  and  Poturalski's  investi¬ 
gations  were  grown  on  cylindrical  surfaces.  The  waves  were 
launched  by  a  fine  wire  located  close  to  the  YIG  surface. 

The  direction  of  propagation  is  noted  by  the  wave  vector  k  in 
Figure  3.1.  These  investigations  were  unique  in  that  no  one 
had  ever  propagated  MSSW's  in  this  fashion.  Conventional 
delay  lines  were  generally  fabricated  from  samples  cut  from 
films  grown  on  1-  to  3-inch  diameter  GGG  flat  wafers. 

Earlier  experiments  in  microwave  magnetics  were 
performed  with  YIG  rods  and  slabs  fabricated  from  bulk 
crystals.  Since  large  defect-free  samples  are  difficult  to 
realize,  a  more  uniform  process  is  desired.  An  outgrowth  of 


the  development  of  magnetic  bubble  memories  was  the  advent  of 
growing  large  homogeneous  thin  films  of  YIG.  The  technique 
used  to  fabricate  these  single  crystal  films  is  liquid  phase 
epitaxy  (LPE).  A  "melt"  is  prepared  with  the  proper  compo¬ 
sition  and  heated  to  just  above  the  saturation  temperature  of 
about  910  C.  When  the  temperature  of  the  melt  is  brought 
below  saturation,  crystallization  occurs.  A  suitable 
substrate  is  brought  to  the  same  temperature  before  being 
immersed  in  the  solution,  and  serves  as  a  seed  for 
crystallization.  The  growth  of  YIG  on  GGG  is  relatively 
fast--on  the  order  of  one  micron  per  minute.  The  GGG  has 
relatively  the  same  lattice  constant  as  YIG.  Pure  YIG  has 
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a  lattice  constant  which  is  about  0.06 %  smaller  than  that  of 
GGG  (Glass  1982:102).  A  defect-free  substrate  generally 
leads  to  good  epitaxial  growth.  Substitutions  may  be  made  in 
the  melt  to  help  compensate  for  the  lattice  mismatch,  but 
these  impurities  also  alter  the  magnetic  properties  of  the 
film.  Therefore,  substitution  is  used  as  a  trade-off  which 
allows  growth  of  thick  films  for  certain  device  applications, 
or  to  provide  temperature  compensation  of  the  material  for 
other  applications  (Glass,  1982:104).  For  example,  it  has 
been  found  that  when  the  lead  content  of  the  melt  equals  the 
platinum  impurity  concentration,  the  melt  produces  films 
which  have  the  narrowest  FMR  (Ferro-Magnet ic  Resonance)  line- 
width  (a  measurement  of  loss  quality)  and  smallest  lattice 
misfit  (Glass,  1976:286). 


SIDE  VIEW 


END  VIEW 


RF  COAX  FEEDS  — ,  RF  COAX  FEED 


\  CENTER  HOLE 
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Figure  3.1  Diagram  of  YIG  Film  Used  for  Rotation 

Rate  Sensing 
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Upon  review  of  the  fabrication  process  of  YIG  films, 
several  reasons  for  anisotropic  characteristics  are  evident. 
One  prospect  is  uniformity.  If  the  film  thickness  is  not 
uniform,  it  has  been  shown  that  magnetic  properties  will 
vary.  Another  reason  for  anisotropy  is  the  stress  induced  by 
the  lattice  mismatch.  Lattice  mismatch  is  the  main  cause  of 
material  defects.  There  is  also  magnetocrystalline 
anisotropy,  which  Poturalski  thought  to  be  the  main  cause  of 
his  poor  propagation  results  (Poturalski,  1979:79).  He  felt 
that  the  cubic  structure  beneath  the  launch  transducer  had  an 
effect  on  the  coupling  efficiency. 


For  normal  delay  line  configurations,  the  relative 
position  of  the  transducer  and  the  film  does  not  change:  but 
with  the  rotation  rate  sensor,  the  change  relative  to  the 
transducer  is  three-dimensional.  This  investigation  will  use 
Poturalski 's  launch  transducers  on  a  flat  sample,  with 
rotation  in  the  plane  of  the  single  crystal  film,  and  thus, 
its  cubic  structure  will  be  available  for  investigation 
(Figure  3.2).  The  Poturalski  transducers  were  chosen  for  two 
reasons.  One,  they  provide  an  ability  to  launch  magnetic 
waves  with  adjustable  coupling,  and  second,  they  provide  high 
isolation  because  of  their  BF  shielding  properties. 


Figure  3.2  Cross-Section  View  of  Poturalski  Probe 
(Poturalski,  1979:42) 

The  experiments  detected  and  recorded  the  influence  of 
the  crystal lographic  structure  on  the  propagation  character¬ 
istics  of  magnetostatic  waves.  These  probes  were  a  key 
factor  in  the  evaluation  of  the  microwave  properties  that 
were  under  investigation.  They  permitted  the  fabrication  of 
a  simple  delay  line  in  the  magnetostatic  surface  wave 
configuration,  a  simple  single  wir?  transducer  to  launch  the 
wave,  and  another  for  the  detection  of  the  surface  wave. 

Normal  evaluation  of  a  delay  line  involves  several 
characteristics,  two  of  which  include  insertion  loss  and  time 
delay.  These  characteristics  will  be  measured  and  recorded 
as  a  function  of  the  rotation  angle  o  (Figure  3.3). 


Figure  3.3  Orientation  of  Disc  to  Rotation 
Angle  6.  (Top  View)  (After  Poturalski,  1979s78) 


As  noted  in  the  previous  chapter,  the  films  under 
investigation  had  (111)  surfaces.  Any  influence  of  the 
crystallographic  structure  should  have  three  high-symmetry 
directions,  as  the  sample  is  rotated  from  0°  to  360°.  There 
should  be  a  sixfold  dependency  present,  and  its  influence  is 
what  this  investigation  is  trying  to  analyze. 

The  samples  to  be  tested  were  provided  by  James 
Sethares  of  the  Rome  Air  Development  Center,  Hanscom  AFB ,  MA. 
The  three  samples  were  grown  by  the  LPE  method.  The  film 
thicknesses  were  11.5,  7.35,  and  7.5  microns  for  the  three 
samples.  Each  had  a  small  lead  impurity,  which  tended  to 
help  the  tensile  misfit?  thus,  the  samples  should  have  very 
little  stress  induced  anisotropy  (Glass,  1982:102). 
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Experimental  Setup 

A  test  apparatus  was  designed  for  accurate  placement 
and  measurement  of  input  and  output  probes.  By  noting  the 
basic  test  configuration  shown  in  Figure  3.3*  it  can  be 
observed  that  the  design  of  the  apparatus  had  to  permit  both 
rotation  of  the  disc  for  the  basic  anisotropic  data,  and 
accurate  movement  of  the  transducers  to  collect  beam  steerin 
data  when  anisotropic  effects  influenced  the  direction  of 
propagation.  The  testing  apparatus  is  pictured  in  Figure 
3.4. 


Figure  3.4  Test  Apparatus 
The  Poturalski  probes  were  attached  to  micrometer 
adjustable  translation  stages.  The  input  probe  varied  only 
in  the  x-direction  in  order  to  optimize  RF  coupling  and 
facilitate  the  insertion  and  removal  of  samples.  The  output 
probe  was  mounted  on  a  three-dimensional  stage.  In  this  way 
beam  steering  could  be  investigated,  as  well  as  time  delay 
and  amplitude  distributions.  The  1-inch  diameter  YIG  films 
with  500  pm  GGG  were  mounted  on  a  i-inch  (6.35mm)  plexiglass 


disc,  and  then  mounted  onto  the  2-inch  high  stand.  The  stand 
was  circular,  with  angles  marked  between  0°  and  360°,  in  53 
increments.  The  machining  of  the  stand  was  accomplished  to 
an  accuracy  of  ±  1  degree  of  arc.  Three-inch  square  perma¬ 
nent  magnets  were  mounted  on  adjustable  slides  in  order  to 
provide  variable  magnetic  bias  fields.  The  permanent  magnets 
provided  a  field  strength  ranging  from  0  to  800  Gauss. 

Pulse  Measurements 

The  circuit  depicted  in  Figure  3.5  was  developed  for 
testing  the  RF  characteristics.  A  continuous  wave  (CW) 
sinusoidal  signal  generator  was  gated  with  a  PIN  diode 
switch.  The  gated  RF  signal  was  split  by  a  10dB  coupler,  and 
the  main  signal  was  passed  through  a  tuner,  for  matching 
purposes,  and  then  to  the  MSSW  delay  line  and  a  TWT  amplifier 
before  detection.  The  coupled  signal  was  used  as  a 
reference . 

From  this  configuration,  data  was  obtained  in  three 
categories,  which  include  signal  power,  time  delay,  and  beam 
steering  of  the  magnetostatic  delay  pulse.  The  YIG  samples 
were  rotated  through  3 6 C  ,  with  data  being  taken  at  10  incre¬ 
ments.  Since  the  single  wire  transducers  produce  a  wideband 
signal,  the  frequency  used  for  testing  was  adjusted  to 
minimize  insertion  loss  when  the  magnetic  bias  field  (H)  was 
fixed  at  300  Gauss.  Evaluation  of  the  results  of  these 
investigations  are  reported  in  Chapter  Four. 
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What  one  observes  from  the  resulting  data  of  the  gated 


RF  measurements  is  some  very  erratic,  if  not  random,  data. 
These  findings  led  to  testing  the  apparatus,  with  the  delay 
line  in  place,  on  a  network  analyzer.  This  procedure 
permitted  a  frequency  domain  picture  of  the  device  to  be 
taken.  Observation  of  the  delay  line's  passband  shifting  up 
and  down  in  frequency,  when  the  film  was  rotated  beneath  the 
transducers,  led  to  testing  the  device  in  an  oscillator 
conf igurat ion . 
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CIRCUIT  DIAGRAM  FOR  PULSE  MEASUREMENTS 

Figure  3.5  Block  Diagram  for  RF  Pulse  Measurements 
Delay  Line  Oscillator 

A  delay  line  oscillator  circuit  was  possible,  since 
there  would  be  no  need  to  change  the  delay  line  configuration 
of  the  test  apparatus  (see  Figure  3-6).  In  this  apparatus. 


Pulse  Measurements 


As  described  in  Chapter  Three,  the  test  fixture  was 
constructed  to  accurately  hold  the  transmitting  and  receiving 
probes  for  profiling  the  characteristics.  Most  of  the 
experiments  used  CW  pulses  at  2400  MHz  (s-band).  Among  the 
first  data  to  be  taken  were  the  field  profiles. 

In  Figure  4.1,  the  three-dimensional  shape  of  the 
amplitude  distribution  can  be  observed.  This  general  shape 
was  true  for  all  the  YIG  samples.  In  Figure  4.2a,  a  cross- 
section  of  the  amplitude  distribution  is  shown.  The  data  of 
Figure  4.2a  was  taken  in  the  YZ  plane  with  z=10mm.  The  hump 
shape  in  the  exponential  falloff  is  due  to  reflection  f-’om 


Figure  4.1  Amplitude  Distribution  of  MSSW  Beam 


the  sample  edge.  Figure  4.2b  shows  the  time  delay  of  a 
5  nsec  wide  pulse  through  the  magnetic  medium  as  a  function 
of  path  length  in  meters. 

From  Figure  4.1,  it  can  be  observed  there  is  a  slight 
steering  of  the  beam  to  the  left.  This  steering  effect 
changed  as  the  sample  was  rotated.  At  times,  the  beam  would 
be  steered  to  the  right,  but  this  did  not  have  a  sixfold 
dependence.  Measurements  of  the  MSSW  beam  in  the  x-y  plane 
provide  a  cross-sectional  view  of  the  time  delay  profile  or 
amplitude  roll  off.  A  set  of  measurements  was  also  taken 
with  Y  =  constant  for  the  beam  steering  profiles. 
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Figure  4. 2. a.)  Amplitude  Profile  b.)  Time  Delay  Profile 


The  resultant  data  was  then  plotted  for  evaluation. 
Some  of  the  resultant  graphs  are  depicted  in  Figures  4.3  and 
4.4.  Data  was  taken  at  10°  increments  and  normalized  to  the 
greatest  value  to  emphasize  any  periodicity.  Visual 
inspection  of  the  data  in  Figures  4.3  and  4.4  shows  some 
random  nature,  and  it  is  very  difficult  to  deduce  anything 
directly.  However,  techniques  do  exist  for  evaluation  of  the 
apparently  random  data.  A  correlation  analysis  was  helpful 
in  interpreting  these  results. 

The  autocorrelation  function  will  be  examined  first. 
The  use  of  this  technique  will  emphasize  any  major  period¬ 
icities  which  may  be  present  in  the  data.  The  autocorre¬ 
lation  function,  R(0),  in  discrete  form  is: 

N-B 

X„X„4S,  e--o»“...3so*  (,,o; 

*  N’6 

In  order  to  validate  the  technique,  a  random  number 
generator  was  used  to  create  test  data  with  a  sixfold  depen¬ 
dency  added  in.  The  data  is  pictured  in  Figure  4.5  with  its 
autocorrelation.  The  correlation  function  was  able  to  iden¬ 
tify  the  sixfold  periodicity.  There  is,  however,  a  drawback 
to  autocorrelation.  This  technique  may  tend  to  hide  one 
periodicity  if  a  more  dominate  periodicity  is  present. 
Therefore,  a  technique  called  cross-correlat ion  is  sometimes 
used.  This  function  is  similar  to  autocorrelation,  but 
instead  of  correlating  data  with  itself,  it  is  correlated 
with  a  known  function  containing  the  dependence  one  is 
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Figure  4.4  Sample  Two.  Pulse  Measurements  of 
a)  amplitude,  b)  time  delay,  and  c)  beam  steering 


searching  for.  Thus,  comparison  of  cross-correlations  will 
reveal  the  relative  strengths  of  the  periodicity. 


One  of  the  main  reasons  for  using  correlation 
techniques  is  evident  from  the  data  in  Figure  4. *4.  The 
amplitude  data  seems  to  have  no  evident  periodicity,  but  when 
the  autocorrelation  is  evaluated,  one  can  see  a  noticeable 
twofold  dependency  present  in  Figure  4.6c.  There  is  also 
some  evidence  of  a  sixfold  dependency.  In  Figure  4.6, 
several  correlation  graphs  are  presented  for  the  experiments 
performed  on  signal  strengths  and  time  delays.  These  graphs 
show  the  general  trend  of  this  data.  They  are  not  very 
convincing,  however,  because  the  sixfold  crystal  symmetry  has 
such  a  large  effect  on  the  magnetostatic  wave. 
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Figure  4.5  a)  Random  Test  Data  with  Sixfold  Dependence  Added 
b)  Autocorrelation 


The  facts  or  conclusions  which  could  be  drawn  from 
these  first  experiments  were  unconvincing  for  a  sixfold 
dependence.  Other  characteristics  were  sought  and  investi¬ 
gated  to  find  if,  indeed,  there  was  any  sixfold  influence. 


Figure  4.7  shows  actual  data  correlation  graphs  of  insertion 
loss,  group  delay,  and  group  velocity.  These  results  were 
also  not  at  all  conclusive.  However,  when  this  data  was 
taken,  the  group  delay  data  was  derived  from  the  phase 
measurement  on  a  network  analyzer.  The  network  analyzer, 
which  could  display  amplitude  and  phase  of  the  transmission 
characteristic  versus  frequency,  revealed  some  interesting 
results.  The  amplitude  versus  frequency  response  illustrated 
the  bandpass  or  "passband"  of  the  delay  line.  When  the  YIG 
sample  of  the  delay  line  was  rotated,  one  could  see  the 
passband  shifting  up  and  down  in  frequency.  This  shifting 
appeared  to  be  related  to  the  internal  DC  field  bias,  because 
holding  0  constant  and  varying  the  DC  bias  caused  the  same 
result.  This  observation  motivated  oscillator  experiments 
which  looked  for  frequency  shifts  related  to  the  crystalline 
an  isotropy. 

Delay  Line  Oscillator  Measurements 

The  circuit  for  frequency  shift  effects  must  be  able  to 
detect  small  frequency  variations  in  the  samples.  One  way  to 
achieve  this  is  to  test  the  sample  as  a  resonator.  The 
optimum  and  easiest  way.  however,  is  to  use  it  as  an 
oscillator.  By  using  the  YIG  sample  as  a  delay  line  in  the 
feedback  loop  of  a  delay  line  oscillator,  the  frequency  of 
oscillation  becomes  sample  dependent.  A  diagram  of  the 
circuit  is  shown  in  Figure  4.8. 
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Figure  4.6  a)  Sample  One,  Time  Delay  Autocorrelation 

b)  Sample  Two,  Amplitude  Autocorrelation 

c)  SampLe  One,  Amplitude  Autocorrelation 
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Figure  M.7  Correlation  graphs  of 

a)  IL-Insertion  Loss 

b)  GD-Group  Delay 

c)  GV-Group  Velocity 


CIRCUIT  DIAGRAM  FOR  DELAY  LINE  OSCILLATOR 


Figure  4.8  Block  Diagram  of  the  Delay  Line  Oscillator 

When  these  tests  were  performed,  everything  was  held 
constant,  except  the  rotation  angle  (o).  The  gain  of  the 
loop  was  adjusted  such  that  only  one  frequency  was 
oscillating.  It  was  easy  to  track  a  particular  wave  number 
in  this  way.  As  the  net  effective  field  in  the  YIG  was 
altered  due  to  the  rotating  sample,  the  frequency  of  oscil¬ 
lation  shifts  slightly.  In  this  way,  measurements  were  taken 
for  F  versus  9,  and  the  results  are  displayed  in  Figure  4.9. 
The  sixfold  dependency  is  easily  identifiable,  even  on  top  of 
a  twofold  dependency.  The  sixfold  dependence  is  more 
noticeable  in  the  second  sample  as  compared  to  the  first 
sample.  The  maximum  frequency  shift  occurred,  however,  in 
the  first  sample  and  it  was  24.35  MHz  (peak  to  valley). 


The  autocorrelations  of  Figure  4.10  show  the  definite 
sixfold  dependency.  The  cross-correlations  of  Figure  4.11 
allow  comparison  of  the  sixfold  dependency  between  samples. 

It  should  be  noted  here  that  the  amount  of  sixfold  dependency 
found  in  each  sample  is  inversely  related  to  the  thickness  of 
the  sample,  but  it  is  not  proportional .  Recall  that  sample 
thicknesses  were  11.5,  7.35,  and  7.5  microns  for  samples  one, 
two,  and  three,  respectively.  Although  this  is  interesting, 
no  solid  conclusion  can  be  drawn  from  it. 

There  was  another  factor  which  also  influenced  the 
testing,  and  it  was  investigated  further.  From  Figure  4.1, 
the  three-dimensional  plot  of  the  amplitude,  it  can  be  seen 
that  the  maximum  amplitude  occurred  slightly  off  the  center 
line  of  the  sample.  It  was  found  that  when  everything  was 
held  constant,  including  the  rotation  angle  (o),  the  location 
of  the  maximum  amplitude  signal  could  be  changed  by  adjusting 
the  CW  frequency  of  the  pulses.  This  beam  steering 
phenomenon  was  recorded  and  examined  in  Chapter  Five. 
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.10  Autocorrelations  for  Figure  M 

a)  Sample  One 

b)  Sample  Two 

c)  Sample  Three 
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.11  Cross  Correlations  of  Figure 

a)  Sample  One 

b)  Sample  Two 

c)  Sample  Three 
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V .  Beam  Steerin 


As  mentioned  in  earlier  chapters,  there  were  experi¬ 
ments  performed  where  the  effects  of  beam  steering  were 
present,  in  addition  to  the  crystallic  anisotropy  influence. 
It  was  found,  for  example,  that  with  all  parameters  held 
constant,  including  0,  that  the  position  of  maximum  signal 
strength  (beam  steering)  could  be  altered  by  simply  changing 
the  CW  frequency  of  the  pulsed  signals.  Some  data  was  taken, 
and  a  plot  of  the  results  is  shown  in  Figure  5.1. 
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Figure  5.1  Beam  Steering  as  a  Function  of  Frequency 
Upon  examination  of  the  test  apparatus,  it  was  found 


that  the  input  probe /transducer  RF  shield  was  not  parallel  to 
either  the  sample  or  the  transducer.  This  was  suspected  to 
be  the  cause  of  the  frequency  dependent  steering  because 


conductors  are  known  to  influence  magnetostatic  waves,  and 
this  non-uniformity  was  probably  causing  some  interference. 
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Figure  5.2  View  of  Non-Parallel  RF  Shield 


To  understand  how  this  physical  perturbation  may  affect 
the  MSW  beam,  we  need  to  view  a  sketch  of  the  geometry. 

Figure  5.3  illustrates  the  regions  close  to  the  transducer. 

The  beam  passes  through  Region  2,  which  represents  the  non¬ 
parallel  shield.  The  velocity  in  this  region  is  a  function 
of  the  spacing  of  the  shield.  Since  the  spacing  is  tilted, 
the  velocity  is  a  function  of  y.  The  frequency  of  the  wave 
remains  the  same,  but  the  velocity  changes,  and  thus,  so  does 
the  wavelength.  A  theoretical  expression  was  sought  to  model 


this  case,  but  none  was  satisfactorily  obtained. 
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The  angle  of  refraction  is  calculated  below: 
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where  Ki  and  K2  are  the  MSSW  wave  numbers  in  their  respective 
regions . 
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Figure  5.4  Diagram  of  the  Oblique  Incidence 
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itself  well  to  probing,  a  semi-circle  of  receivers  was 
arranged  approximately  1  inch  from  the  input.  This  apparatus 
provided  for  a  more  stable  and  controlled  experiment. 


Figure  5.5  Test  Apparatus  for  Examining  Beam  Steering 

Since  magnetostatic  waves  are  dispersive,  it  is  easily 
confirmed  that  it  is  the  velocity  differential  across  the 
boundary  that  can  steer  the  MSW  beam.  The  dispersive  nature 
propagates  different  frequencies  at  different  velocities; 
therefore,  different  frequencies  will  be  steered  by  differing 
amounts.  Figure  5.6  demonstrates  this  property.  Two  pulses 
of  different  CW  frequencies  were  propagated  through  the 
triangular  metallization  region.  While  the  pulse  of  a  3-  10 
GHz  signal  was  steered  to  receiver  number  four,  the  pulse  of 
a  3.55  GHz  signal  was  received  at  channel  three.  The  first 
demonstration  of  the  property  of  refraction  in  magnetostatic 
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waves  was  performed  with  this  device.  Thus,  it  has  been 
shown  that  the  metallization  does  indeed  have  an  affect  on 
the  magnetostatic  wave.  This  is  convincing  evidence  which 
supports  the  conclusion  that  the  tilting  RF  shield  did 
influence  and  bias  the  measurements  of  the  probe  apparatus. 


SEPARATION  OF  TWO  SIMULTANEOUS 
CW  PULSES 


Figure  5.6  Separation  of  Pulse  with  Differing  CW  Frequencies 


Conclusions  and  Recommendations 


This  thesis  investigated  anisotropic  effects  in 
magnetostatic  surface  waves.  There  were  several  experiments 
performed  during  this  investigation,  and  important  conclu¬ 
sions  may  be  drawn  from  them.  The  first  experiments  were 
simple  delay  line  measurements  of  samples  at  variable  input 
and  output  transducer  locations.  The  samples  were  rotated 
beneath  the  transducers,  while  data  for  signal  attenuation, 
time  delay,  and  beam  steering  were  collected.  These  measure¬ 
ments  resulted  in  the  collection  of  what  appeared  to  be  very 
random  data.  Auto-  and  cross-correlation  techniques  were 
used  to  analyze  the  data.  It  was  found  that  there  were 
indeed  periodicities  present  in  the  data.  A  twofold 
periodicity,  or  every  1 80  ,  was  observed  in  some  of  the  data. 
This  was  expected,  since  whenever  the  sample  was  mounted  at  a 
slight  tilt  from  the  horizontal,  the  twofold  dependency 
resulted  from  the  fact  that  when  rotated  1 8 0 ",  the  sample 
would  be  at  the  same  tilt  angle,  but  propagation  would  be  in 
the  opposite  direction.  Even  with  random  data  and  a  dominant 
twofold  dependency,  a  sixfold  crysta 1 1 ograph i c  dependency  was 
detected . 

The  second  set  of  experiments  involved  measurement  of 
the  YIG  samples  using  a  network  analyzer.  For  a  delay  line 
configuration,  insertion  loss,  group  delay  and  group  velocity 
data  were  recorded.  All  of  the  data  was  found  to  be 
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consistent  with  known  MSW  characteristics.  It  was  from  these 
measurements  that  a  frequency  shifting  effect  was  noticed. 

The  delay  line  oscillator  measurements  revealed  the 
most  informative  crystallographic  results.  This  config¬ 
uration  was  more  sensitive  to  the  anisotropic  effects,  which 
change  the  effective  H  field  acting  on  the  magnetic  spins  of 
the  thin  ferrite  film.  The  frequency  of  oscillation  is 
directly  dependent  on  the  applied  and  effective  H  field.  The 
oscillation  frequency  was  shifted  slightly  by  a  relative 
rotation  between  the  sample  and  the  transducers.  This  shift 
was  shown  to  have  sixfold  symmetry  through  a  rotation  of  360  , 
The  amount  of  shift  is  in  close  agreement  with  anisotropic 
constants  of  YIG  measured  by  other  methods.  A  theoretical 
shift  for  if,  using  these  anisotropic  constants,  was 
calculated  to  be  19  MHz  as  compared  to  the  measured  results 
of  24  MHz  shown  in  Chapter  Four.  These  results  represent  a 
third  technique  for  the  measurement  of  the  magnetocrystalline 
anisotropy  and  show  agreement  with  the  results  of  Sethares 
and  Tsai  (Sethares,  1977:1236),  and  their  ferromagnetic 
resonance  technique  and  Hanna,  et  al,  for  their  surface 
acoustic  wave  phase  shift  technique  (Hanna.  1983:1803).  The 
oscillator  experiment  is  the  first  positive  verification  of 
anisotropic  effects  on  the  propagation  of  magnetostatic 
surface  waves. 

The  work,  which  studied  the  beam  steering  effects  of 
the  Poturalski  probes,  helped  produce  the  most  rewarding 
results.  As  a  result  of  the  findings,  a  second  test  fixture 


was  constructed  to  verify  the  amount  of  influence  of  a 
closely  spaced  metal  conductor.  The  geometry  of  the  fixture 
enabled  the  first  demonstration  of  magnetostatic  wave 
refraction . 

In  summary,  this  investigation  verified  the  existence 
of  anisotropic  effects.  It  was  shown  that  the  amount  of 
influence  is  negligible  in  delay  lines.  However,  for 
applications  which  are  highly  freqjency  dependent,  such  as  an 
oscillator  or  a  phase  shifter,  there  is  a  need  to  include 
magnetocrystalline  constants.  The  results  of  this  thesis 
agree  with  measured  magnetocrystalline  constants  of  Everett 
and  Hughes  (Everett,  1970:1120).  The  delay  line  oscillator 
results  are  a  new  way  to  identify  and  predict  the  influence 
of  magnetocrystalline  anisotropy  on  magnetostatic  surface 
waves . 

Re commend at  ions 

A  recommendation  for  further  study  concerning  aniso¬ 
tropic  effects,  would  be  to  devise  a  more  accurate  scheme  for 
mounting  the  sample  transducers,  which  would  make  their 
contributions  of  mounting  error  less  significant.  This 
would,  hopefully,  reduce  some  of  the  random  fluctuations  of 
data  taken  for  the  delay  line  signal  strengths,  time  delays, 
and  beam  steering.  Another  suggestion  would  be  to  use  a 
narrowband  transducer  to  generate  fewer  frequencies  for 
propagation  across  the  sample.  This  would  enable  the 
experimenter  to  more  accurately  track  the  signal  path  and  its 


characteristics . 


The  work,  which  explored  the  phenomena  of  MSW  refrac¬ 
tion,  shows  promise  as  a  frequency  separator  and  should  be 
investigated  further.  It  should  be  possible  to  exploit  the 
features  of  refraction  for  frequency  selective  and/or 
dependent  devices  for  use  in  electronic  warfare  counter¬ 
measure  applications.  A  U.S.  patent  application  for  such  a 
device  has  been  filed  (Floyd,  1984). 
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